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I argue that in contrast to B < B^ , where the Vortex Liquid ( VL) freezes into the Strongly-pinned 
Bose Glass phase (SBG), with flux-line vortices localized by the columnar defects, for B > B,p, the 
additional vortices see a significantly weaker random ^-independent potential due to the caging by 
the vortices strongly pinned on the columnar defects. The result is that for B > B$, upon cooling, 
the VL undergoes a sharp crossover to the Interstitial Liquid (IL) where SBG coexists with a resistive 
liquid of interstitial vortices. At lower T the IL freezes into a distinct Weakly-pinned Bose Glass 
phase (WBG). Therefore for B > B^ the true linear superconductor only forms at T well below the 
B < B$ irreversability line. 
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In high-T c superconductors, because of the enhanced 
fluctuations, as compared to conventional superconduc- 
tors, over a large portion of the B-T phase diagram, the 
Abrikosov lattice melts [1], [2] into a resistive Vortex Liq- 
uid (VL) of highly mobile flux lines. At lower tempera- 
ture and magnetic field, below the irreversibility line, de- 
pending on the nature of the disorder, the VL freezes into 
one of the proposed new glassy phases, leading to true 
linear superconductivity. In the presence of point dis- 
order such as oxygen vacancies and interstitials, Vortex 
Glass (VG) [3] has been proposed as the low-temperature 
phase replacing the Abrikosov lattice. 

With the success of Civale, et al. [4] in significantly 
shifting the irreversibility line to higher T and B by ir- 
radiating the superconductor with heavy energetic ions, 
much of recent attention has focussed on flux-line vor- 
tices in the presence of correlated disorder such as colum- 
nar defects. It was previously known that the inter- 
acting flux lines map onto the quantum mechanics of 
two-dimensional bosons [1]. Key progress was made 
when Nelson and Vinokur (NV) realized that the random 
columnar defects naturally map onto a time-independent 
random potential for the bosons [5], resulting in the 
previously studied problem of two-dimensional quantum 
bosons on a random substrate [6]. NV interpreted the 
irreversibility line as the phase transition where the VL 
freezes into a superconducting Bose Glass (BG) of vor- 
tices localized on columnar defects. The mapping led to 
a detailed analysis of the BG transition and of the prop- 
erties of the low temperature BG phase. 

Most of the analysis of Ref. [5] had focussed on the 
regime B < B^ with the matching field B^ = 4> /d 2 cor- 
responding to an average spacing d between the colum- 
nar defects. In their work, however, NV did point out, 
that since the random potential for B > B$ is also z- 
independent, the BG will exist in this regime as well. In 
this Letter, I will reexamine in more detail this local- 
ization problem, contrasting the regimes B < B^ and 
B > B$. I will argue for a quantitatively modified phase 



diagram of the BG superconductor for B > B^, with a 
potentially sharp crossover as a function of T or a 1st- 
order transition, from the VL phase to the Interstitial 
Liquid (IL), and a strong suppression of the BG transi- 
tion (the irreversibility line) to lower T, down to where 
the melting line (in the absence of columnar defects) used 
to be (see Fig.l). Although I do not offer many new tech- 
nical details (most of these can be found in Refs [6,5]), my 
arguments have concrete physical consequences that are 
experimentally verifiable. The important consequence of 
my work is that the introduction of columnar defects 
should be much less effective in increasing the pinning, 
raising the critical currents, and shifting the irreversibil- 
ity line to higher B and T in the B > B^ regime, in 
contrast to the B <B$ regime. 

The simplest standard model that captures the physics 
of N interacting flux-line vortices through a sample of 
thickness L, in the presence of columnar defects, is de- 
scribed by the free energy F, 
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where r, (z) is the local deviation of the «th vortex from 
the applied field ^-direction, V(r) = 2eoKo(r /\ a b) is the 
interaction, which is logarithmic up to the London pene- 
tration length X a b, and falls off exponentially for r > A„6. 
U(r) is the random potential of range b « 70A, induced 
by the columnar defects, with disorder-averaged vari- 
ance U(f)U(?) = 2U^{b A /d 2 )5^\r-r l ). The line ten- 
sion e; = (M a i,/M z )eo ln(A &/£) is significantly reduced 
relative to the condensation energy eo = ((/>o/47rA a (,) 2 
due to the considerable anisotropy of high-T c materials 
(M a b/M z « 1/100). 

As discussed in detail by NV [5], building on the ideas 
of Ref. [6], the above model for B < B^ predicts three 
distinct vortex phases. I first remind the reader of the 
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properties of these phases and then turn to the new pre- 
dictions, relying on many results of NV. 

At high temperature and magnetic field, the fluctu- 
ating vortices can move freely through a sample, "dif- 
fusing" along the z-axis, < \fi(z) — ?i(z)\ 2 > « 2zT/ei. 
This leads to a resistive VL phase (corresponding to the 
superfluid phase of the quantum bosons) , with the flux- 
flow resistivity described by the Bardeen-Stephen (BS) 
formula, Pff = p n B/B C 2; p n is the normal-state resistiv- 
ity. The bulk modulus en and the tilt modulus C44 ~ eiB 
are finite, while the shear modulus vanishes in this phase. 
For B < B$, as T or B is lowered below T BG (B), the BG 
transition takes place, at which the flux lines localize on 
the columnar defects, as described by the finite transverse 
wandering length l± and < \fi(z) — fi(z)\ 2 > = l\. Upon 
approaching the transition the tilt modulus (proportional 
to the inverse superfluid density for the quantum bosons) 
will diverge leading to the transverse Meissner effect, re- 
quiring a finite tipping angle 6 C ss (Uob 2 n/ 'dPc^) 1 / 2 of 
the external field (relative to the columnar defects) to 
change the average internal magnetic field direction. 

At the BG transition, postulating two divergent length 
scales, the localization length l± ~ 1/(Tbg — T) v , the 
correlation length along the ^-direction 2y ~ l\, and the 
corresponding divergent relaxation time r ~ l\_ , one can 
construct the BG scaling theory. This analysis predicts 
the I-V curve given by E(J) = ZJ 1 " z 'f +/ _(JZ_ l Z||0 o /cT), 
with the scaling function F + /_ (x), which in the VL phase 
(T > Tbg) is F + (x) ~ x, leading to a finite linear resis- 
tivity pff ~ (T — TbgY ^ z ~ 2 ^- I n the truly supercon- 
ducting BG phase, p// vanishes and various transport 
regimes lead to highly nonlinear I-V characteristics. At 
small J, in the variable-range-hopping regime, where the 
vortices hop to distant columnar defects via superkinks, 
E(J) ~ exp(— 1/ J 1 / 3 ). At high J, closer to the critical 
current J c , E(J) ~ exp(— 1/J), and at the BG transition 
E(J) ~ '/ 3 . The scaling theory also predicts the 

existence of 9 C ~ (Tbg — T) 3 " , which leads to a cusp 
in resistivity R(0) and in the Tbg(H±), as is observed 
in experiments on crystals with correlated pinners [9,10], 
the features that distinguish BG from the isotropic VG 
phase. Finally, as pointed out by NV, a highly correlated 
Mott Insulator (MI) phase might appear asB4 B^, in 
which both C44 and C\\ are infinite. 

I now closely examine the situation for B > B^ (see 
Fig.l). Because of the divergence of the bulk modu- 
lus en, the MI phase is expected to exhibit rigidity to 
the introduction of any additional flux lines, and in this 
way behaves similarly to the Meissner phase. In fact, 
recent magnetization-relaxation experiments appear to 
verify its existence, observing a drastic suppression of 
the magnetization relaxation rate for B « B$ [8]. 
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FIG. 1. The proposed phase diagram for the BG supercon- 
ductor. As discussed in the text, the SBG phase, where the 
vortices are strongly pinned by columnar defects, is replaced 
by two phases: the resistive IL and the superconducting WBG 
phase for B > B^. The superconducting MI phase is the 
Meissner phase for the interstitial vortices. 

As the applied field is increased above B$, additional 
vortices begin to enter the sample. It is energetically 
more favorable for these vortices to occupy the intersti- 
tials between the vortices already strongly pinned by the 
columnar defects. At very low temperatures these inter- 
stitial vortices will also be localized by the repulsion V(r) 
from the vortices on columnar defects, caged-in by the 
nearest neighbors. Since the randomly localized colum- 
nar vortices will appear to the interstitial flux lines as a 
random, ^-independent potential, I expect that the in- 
terstitial vortices will also freeze into the BG phase [5]. 

The resulting Weakly-pinned Bose Glass (WBG) phase 
will have a vanishing linear resistivity, and the transport 
describable by the BG scaling theory. In fact, the equi- 
librium properties of this phase should be similar to that 
of the SBG, described above, in all respects except for 
the pinning energy scale which determines Tbg(B) and 
J c . In the SBG phase (at low T « T^ BG ss 70 K, where 
the effects of thermal renormalization are not impor- 
tant), the columnar defect pinning energy is determined 
by C^o ~ eo with the corresponding J^ BG = cUo/(fiob. 
Quantitatively I expect that the interstitial vortices con- 
trolling the formation of the WBG phase are significantly 
weaker localized. Because the repulsive energy V(r) and 
the columnar defect localization energy U (r) are both of 
the order of the condensation energy eo, it is difficult to 
have a good model-independent estimate of the localiza- 
tion energy scale of the interstitial vortices. One possible 
guess is Ui ~ li/dn 1 / 2 « Uq, which leads to a reduced 
jWBG = tfji/fod « Jf BG . Physically it is clear that 
Ui « Uq because of the competition between the repul- 
sion from the neighboring vortices strongly localized on 
columnar defects and the "attraction" by the columnar 
defects. One limit in which it is unambiguously clear that 
Ui « Uq is when the spacing between the columnar de- 
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fects d » X ab , in which case £/, ss U e~ d l Xah «U Q . In 
fact, recent experiments of Ref. [8], which infer from the 
size of the hysteris loop two distinct critical currents for 
B < B and B > B^ with jf BG « 2J C W/BG , give strong 
experimental support to my arguments. The existence 
of two critical currents will also lead to a double-sloped 
spatial magnetic field profile of the Bean critical state [7] . 

Given this separation of the pinning energy scales for 
B > B^, in exploring the physics beyond linear response 
it is convenient to treat the vortex system in a two-fluid 
model spirit, with the WBG phase consisting of vortices 
strongly pinned on the columnar defects and a fraction 
of interstitial vortices localized by the repulsive interac- 
tions. This two-fluid model then leads to the shape of the 
I-V curve which is significantly different from the SBG 
phase for B < B^. As illustrated in Fig. 2, the I-V curve 
is characterized by two depinning transitions at J^ BG 
and Jf BG to the flux-flow regime approaching that of 
the IL and VL, respectively. 




FIG. 2. The expected I-V curve in WBG phase, illustrating 
two depinning transitions, characterized by J™ BG and J^ Ba 
and corresponding asymptotic flux flow resistivities. The in- 
set shows qualitatively J C {B), observed in Ref.8. 

The above arguments imply that for B > B$, upon in- 
creasing temperature, the weakly-pinned interstitial vor- 
tices will be first to delocalize. When their transverse 
wandering exceeds their separation the WBG phase will 
melt (most likely via the BG transition) into the resistive 
IL. The resulting IL phase will have properties interme- 
diate between that of the SBG and the VL, calculable 
within the simple two-fluid model. In this regime the in- 
terstitial vortex fraction is completely delocalized, while 
the B^/B fraction of flux lines remains pinned along the 
columnar pinners for a significant longer distance along 
the defect. This distinction between the strongly pinned 
columnar vortices and the interstitial liquid vortices is 
blurred by the derealization of columnar vortices at long 
lengths and time scales and by exchange processes be- 
tween the vortices. Therefore from the point of view of 
the symmetry of the order parameter the IL is identical 



to the VL phase and two are just different regimes of the 
same normal metallic phase. Another real possibility is 
that a 1st order transition separates IL phase from the 
VL phase, which would lead to a jump in quantities (e.g. 
linear resistivity) that depend on the density of mobile 
vortices. In either case, the IL will exhibit quantitatively 
observable differences in behavior from the usual VL. The 
linear flux-flow resistivity in the IL is determined primar- 
ily by the interstitial vortex fraction, and is given by the 
generalized Bardeen-Stephen formula 

pf s = Pn (B - B^)/B c2 (2) 

It is possible that fluctuations will modify the dependence 
on B — B$ to another power law, but pj L will vanish as 
the MI phase is approached. 

As I discussed above, in the IL a mobile fraction 
(B — B$)/B of vortices is in the interstitials and the 
remaining vortices B^/B are localized on defects, with a 
constant dynamic exchange between the two "species". 
Since the vortices are really identical, each flux line 
spends the corresponding fractions of time in the inter- 
stitials and localized on columnar defects. During the 
fraction of time that the flux line is localized on a de- 
fect no diffusion takes place. Therefore the transverse 
"diffusion" coefficient D z = T/e;, determining the trans- 
verse wandering length l±(z) = \/2D z z, is significantly 
renormalized in the IL regime to D\ L = D Z (B — B^)jB. 
Given the fact that the tilt modulus e/ maps onto the in- 
verse of the superfluid density in the corresponding prob- 
lem of quantum bosons I find that in the IL the corre- 
sponding superfluid density appears to be reduced by the 
(B — B^)jB fraction. The phenomenology of the IL is 
therefore quite similar to the superfluid Helium on a ran- 
dom substrate. There one observes that the low T regime 
of the superfluid phase appears to have a substantially 
reduced superfluid density as compared to the extrapo- 
lation from the high T regime and the observations are 
well described in terms of a superfluid fraction of He 4 on 
top of a "dead layer" of nonsuperfluid He 4 . The mapping 
of the IL onto the 2d quantum bosons therefore maps the 
interstitial vortices onto the superfluid fraction and inter- 
prets the observed "dead layer" in terms of the vortices 
strongly pinned on columnar defects. I therefore expect 
that this phenomenon should be observable directly in 
the flux-line vortex system. 

It is expected that Bose Glass can be melted by tip- 
ping the applied field relative to the columnar defect 
^-direction [5]. However, unlike the case for B < B^, 
two critical tipping angles Q^bg ^ ^Jy^Q _ ~B$) « 
9f BG oc \/UoB will appear, corresponding to Twbg(H±) 
and Tsbg(H±) phase boundaries. At the lower boundary 
the WBG melts into the IL, with the interstitial vortex 
fraction v = (B — B^)fB depinning first, followed by the 
complete depinning of the columnar vortices at Osbg, 
crossing over (or making a 1st order transition) into the 
VL. In the IL phase I expect to see striking features in 
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the resistivity as a function of tipping angle, p(0), for J 
perpendicular to z (similar signatures will appear in the 
WBG phase). For 8 < 8f BG , p{8) will follow the usual 
p(8) « vpjj j sin 2 (8) behavior, which will abruptly switch 
to a higher curve p(8) ss pj^sin 2 (8) for 6 > 8f BG , as il- 
lustrated in Fig. 3. 
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FIG. 3. The expected dependence of resistivity in IL phase, 
as a function of the tipping angle with J ± z (full curve). The 
dot-dashed curve is the smooth sin 2 (8) dependence of the IL 
fraction and the dotted curve is the response of the vortex 
fraction which is pinned by columnar defects up to the criti- 
cal angle 6»f SG . 

Unlike the transition for B < B$, for B > B$, 
starting from the VL and lowering temperature, I pre- 
dict a potentially sharp crossover or 1st order transi- 
tion to the IL (at Til), followed by the BG transition 
to the WBG (at Twbg)- Given the two pinning energy 
scales, the transition temperatures can be roughly esti- 
mated as the temperature at which the thermally renor- 
malized localization length l±(T) ss Iq(T/T*) 2 reaches 
the order of the corresponding vortex line separation 
n -1 / 2 , where lo is d and b for interstitial and colum- 
nar vortices, respectively. The thermal crossover tem- 
peratures in the two phases for B > B^ are given by 
Twbg « (2e,ePU i ) 1 ^ and T* SBG « (2e,6 2 [/ ) 1/2 leading 
to T WBG (B) « T^ BG (0 o /d 2 (B - B^)) 1 / 4 and T IL (B) « 
TsBoi&o/tfB) 1 / 4 . The transition from the VL to the IL 
will manifest itself as an additional plateau in the resis- 
tivity p(T) which will disappear for B < B$. 

It is clear that my arguments strongly rely on the sep- 
aration between the interaction [/, and columnar Uq pin- 
ning energy scales. I therefore present additional exper- 
imental evidence supporting the physical picture advo- 
cated above. Before the irradiation by heavy ions (i.e. 
in the virgin sample without columnar defects) the ir- 
reversibility line was the freezing/melting line into the 
Abrikosov lattice (possibly modified by weak point dis- 
order at very long scales, which I ignore here). This freez- 



ing transition is determined by the vortex interaction en- 
ergy scale which I expect to be similar to However, 
the experiments of Civale, et al. [4] find that the melt- 
ing line is significantly shifted to higher temperature and 
magnetic field after the introduction of columnar defects, 
with the new energy scale determined by the columnar 
defect pinning energy Uq. I believe that this observed 
shift is the lower bound on the difference in pinning en- 
ergy of the interstitial and the columnar vortices. In this 
sense, for B > B$, the Twbg(B) transition from the 
IL to the WBG is the resurrection of the virgin melting 
line, smeared to a second-order transition by the colum- 
nar disorder. 

In summary, I have presented arguments for the ex- 
istence of the Weakly-pinned Bose Glass and the Inter- 
stitial Liquid phases for B > B$, and described their 
phenomenology, which is quantitatively quite different 
from the Strongly-pinned Bose Glass and the conven- 
tional Vortex Liquid phases in the B < B^ regime. I 
have discussed some of the experimentally observable fea- 
tures in the IV-curves and in tilting experiments that are 
unique to these new regimes. The observation of a shift 
of the irreversibility line to higher T and B for B < B^ 
when columnar defects are introduced, and the lack of the 
corresponding shift for B > B§ is a direct consequence 
of the presented picture. 
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